A multiphysics and multiscale software environment for modeling astrophysical systems

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly-coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a "Noah's Ark" milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multi-scale and multi-physics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.liComment: 24 pages, To appear in New Astronomy Source code available at http://muse.l

Stellar dynamics in young clusters: the formation of massive runaways and very massive runaway mergers

In the present paper we combine an N-body code that simulates the dynamics of young dense stellar systems with a massive star evolution handler that accounts in a realistic way for the effects of stellar wind mass loss. We discuss two topics: 1. The formation and the evolution of very massive stars (with a mass >120 Mo) is followed in detail. These very massive stars are formed in the cluster core as a consequence of the successive (physical) collison of 10-20 most massive stars of the cluster (the process is known as runaway merging). The further evolution is governed by stellar wind mass loss during core hydrogen burning and during core helium burning (the WR phase of very massive stars). Our simulations reveal that as a consequence of runaway merging in clusters with solar and supersolar values, massive black holes can be formed but with a maximum mass of 70 Mo. In small metallicity clusters however, it cannot be excluded that the runaway merging process is responsible for pair instability supernovae or for the formation of intermediate mass black holes with a mass of several 100 Mo. 2. Massive runaways can be formed via the supernova explosion of one of the components in a binary (the Blaauw scenario) or via dynamical interaction of a single star and a binary or between two binaries in a star cluster. We explore the possibility that the most massive runaways (e.g., zeta Pup, lambda Cep, BD+433654) are the product of the collision and merger of 2 or 3 massive stars.Comment: Updated and final versio